Bicocca-FT-03-7 
DCPT/03/32 
DESY 03-027 
IPPP/03/16 



o , 

O; On a possible measurement of as from B-B 

^: CORRELATIONS IN DECAY 



X 



A. Brandenburg^, P. Nason^ and C. Oleari^ 



<: 

(N 

^ DESY-Theorie, 22603 Hamburg, Germany 
> ; "^INFN, Sezione di Milano, Piazza della Scienza 3, 20126 Milan, Italy 

^ ; ^Department of Physics, IPPP, South Road, Durham DHl 3LE, UK 

(N 

O 
m 
o 



Abstract 



Oh! 

r~| , Motivated by recent preliminary results from the SLD Collaboration on the 

measurement of angle-dependent B-B energy correlations in — > bb events, we 
propose a class of observables that can be computed as a power expansion in the 
strong coupling constant as , of order as at the Born level and that can be used for a 
I precision measurement of as(M^). We compute their next-to-leading order O (af) 

corrections in the strong coupling constant, including exactly quark-mass effects. 
We show that, in the theoretical evaluation of these quantities, large logarithms of 
the ratio of the mass of the final quark over the centre-of-mass energy cancel out. 
Thus, these variables have a well-behaved perturbative expansion in as{Mz)- We 
study the theoretical uncertainties due to the renormalization-scale dependence and 
the quark-mass scheme and we address the question of which mass scheme is more 
appropriate for these variables. 



1 Introduction 

Hadronic final states in e^e~ anniliilation have been intensively studied in recent 
years, both at LEP and SLC, in order to perform QCD tests and to determine as [T]. 
Interesting studies on bottom flavoured hadronic events have also been performed, in 
order to test the flavour independence of the strong coupling constant and to search for 
"running mass" effects |2l El IH El El EI • However, up to now, heavy flavoured jets have 
been studied with the same tools used to study light flavoured jets. For example, the same 
shape variables used for light-quark jets (like thrust, cluster multiplicities, etc.) have been 
applied to heavy flavoured events. On the other hand, heavy flavoured events may have 
some advantages, because the kinematic distributions of the produced heavy quark can 
be computed theoretically with much better accuracy than the corresponding light-quark 
quantities. Furthermore, the kinematic of the final-state heavy flavoured hadron is more 
closely related to that of the final-state quark, since hadronization effects involve energy 
scales that are much smaller than the quark mass. 

Recently, the SLD Collaboration [8 has published preliminary data on the double 
inclusive B-B production, and studied quantities (proposed in Ref. ^9J) that carry infor- 
mation about the angular dependence of the B-B energy correlations. 

The main purpose of the present work is to demonstrate that a slight variant of the 
quantities introduced in Ref. P has suitable properties to perform a novel, heavy-flavour- 
specific determination of the strong coupling constant. We will demonstrate the following 
facts: 

• The observables we are considering have a well-defined, finite expansion in the strong 
coupling constant. 

• Non-perturbative corrections of order Aqcd/"^, where m is the heavy-quark mass, 
cancel in these variables. The only remaining non-perturbative corrections are of 
order Aqcd/£'cm, where i?cM is the centre-of-mass energy. 

• Potentially large terms, enhanced by powers of log(£'cM/''^), cancel at all orders in 
perturbation theory, provided ag is evaluated at a scale near -Ecm- This cancella- 
tion holds as long as one can neglect multiple heavy-flavour pair production in our 
process, which is certainly the case for h production in decays. 

The quantities we are considering have a well defined massless limit, provided one can 
neglect multiple heavy-flavour pair production. This limitation is what makes these quan- 
tities heavy-flavour specific. In fact, only for heavy flavours multiple pair production is 
truly negligible. Thus, although these quantities have formally a well defined massless 
limit, they cannot be used in a straightforward way to study light-flavour production. 

Tools to compute heavy-flavour production at the next-to-leading-order (NLO) level 
have been available in the hterature for the past few years [TUl HH [121 UHl El- We 
have used our calculations O CSl Ull in order to compute the NLO corrections to the 
observables we propose. 



1 



In Sec. we define the variables we are considering and show that their perturbative 
expansion in terms of the strong coupling constant as (Mz) is well behaved, and free from 
large logarithms of the ratio of the quark mass over the centre-of-mass energy. The SLD 
Collaboration has provided us with (preliminary) data for these new variables and 
thus we have been able to check our NLO predictions against the data. 

In Sec. 121 we investigate the massless limit of these variables, and show explicitly that 
potentially large logarithms cancel out. 

In Sec. Hlwe study theoretical uncertainties of our results. We analyze the dependence 
of the NLO results on the renormalization scale and compare the results obtained with 
the pole-mass definition to those obtained using the MS mass, addressing the question of 
which of the two mass definitions is more suitable for this calculation. Finally, in Sec. El 
we give our conclusions. 

In Appendix we collect some leading order (LO) formulae used in our analysis, and 
in Appendix El we comment on the renormalization scheme used in this paper. 



2 Double-inclusive heavy-quark cross section 

Heavy flavoured hadron production proceeds through the process 

e+e-^Zy-f^Q + Q + X , 

(where Q is a heavy quark of mass m) followed by the hadronization of the final heavy 
quarks. In the following we will focus upon B meson production. 

Defining the B hadron scaled energy Xg = 2Eg/Ec^, we consider the following quan- 
tities ^ 

A = - fx^~'^dxg, (2.1) 



cr 



B 



G,W ^ (2,3) 

where dag/dxE is the differential cross section for the final-state hadron production, is 
the angle between the two hadrons, ag the total cross section and Xg, ,„ are the scaled 

1/2 

energies of the two final B hadrons. 

The quantities studied by the authors of are related to the ratios Gij/Gu by a 
calculable perturbative factor {PiPj/Pf in the notation of [HI). In Ref. [8^ these ratios 
were compared to experimental results. However, these quantities have a few drawbacks. 
First of all, they depend upon a "theoretical" factor. It would be preferable to deal with 
quantities that have an experimental definition, free of theoretical assumptions. Second 
and most important, they have a perturbative expansion that starts at leading order with 

"'^The quantities Di and Dij were also considered in 9 and the ratios Gy in |16| . 
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a constant. Therefore, even when corrected at the NLO level, they would allow only for a 
LO determination of as- Conversely, they can be used to test the production mechanism, 
since in the ratios as cancels at first approximation. 

In this work, we claim that the quantities Gtj, defined in Eq. ()2.3|) . should instead be 
studied in order to perform precision QCD tests. We will in fact show in the following that 
these quantities have a perturbative expansion in q;s(-E'cm) with coefficients that remain 
finite even in the limit of large E^yi/m ratios. 

We write the hadronic differential cross sections as 

= [ ^'D^p{x)5{xb- xxb)dxbdx, (2.4) 

dXg^dxB^d COS (p J dXb^ctXfcj" cos 

X 5{xb^ — xiXh^) 5{xb^ — X2Xb2) dxb^dxh2dxidx2 , (2.5) 

where dcxb, Xh, Xh^^^ are the quark differential cross section and the quark scaled energies. 
The quark differential cross sections dcxb are calculable order-by-order in perturbative 
QCD, since the heavy-quark mass m acts as a cut-off for final-state collinear singularities. 
However, in order to go from the QCD partonic cross section dab/dxb to the hadronic one 
dag/dxg, we have to take into consideration non-perturbative effects, of order Aqcd/""^, 
that are present in the hadronic cross section. We assume that all these effects are 
described by a non-perturbative fragmentation function P^p- 

The quark differential cross sections, although calculable, contain logarithmically en- 
hanced terms of the form log' {m/ E^y^, so that their perturbative expansion is not well 
behaved for Ecm 3> m. 

In the quantities Gij{(j)) these large logarithms cancel, together with the dependence 
upon the non-perturbative fragmentation function. This is a consequence of the factor- 
ization theorem, if one makes the additional assumption that secondary production of 
heavy-flavour pairs is negligible, which is in fact the case at the Z° peak. 

The factorization theorem states that mass singularities in inclusive cross sections can 
be absorbed into process independent, universal factors. In our case, mass singularities 
are precisely the log^Ecm/fn) terms. Thus, for Ecm m we have 

= [ ^Vb{z) 6{xb - zxb) dxbdz , (2.6) 
dxb J dxb 

;Vb[Zi)Vb[Z2) 



dxb^dxb2d COS (j) J dxbj^dxb^d cos I 

X 6{xbi — ziXbJ 6{xb2 — Z2Xb2) dxb^dxb2dzidz2 , (2.7) 

where the hatted cross sections do not depend upon m, and all terms that are large in 
the m— >-0 limit (i.e. terms proportional to powers of log (-Ecm /"^)) are absorbed in the 
fragmentation functions Vb{z) ^. Inserting Eqs. (|2.6|) and (j2.7p into Eqs. (|2.4p and (j2.5|) 

^Notice that if we had not assumed a negligible secondary heavy-flavour pair production, wc should 
have introduced also the gluon and light-quark component of the heavy-quark fragmentation function. 
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and considering the definitions in Eqs. ()2.H) and ()2.2|1 . we get 



X 



T>^,Jx) dx 



'VJz)dz] — I xl-^^dx 



D^M) 



X 



x\ ^V^p{xi) dxi 
J zl'Vtizi) dzi 



i-i 
X2 



^P^X2) dX2 



dxb 



{21 



V 

Z2 -Vb{z2)dz2 



X ■ 



0-f 



b2 



dxh^dxh^d cos (p 



(2.9) 



and we see that, in the ratio defining Gij{(f)) (see Eq. ()2.3|) ). the factors (that contain 
mass singularities) cancel out, leaving a finite expression ^. 

Since non-perturbative effects cancel in the definition of Gij, from now on we will make 
no distinction between the quantities Di, Dij defined at the quark and at the hadron level. 
With the available codes for the calculation of NLO corrections in the production of heavy 
quarks [131 Ej at e~^e~ colliders we can compute the contributions up to order to 
the quantities Di, Dij{(j)) and Gij{(j)). We can write the perturbative expansions for Di 
and Dij{(j)) for massive final-state quarks, renormalized in the pole-mass scheme, in the 
region away from the two-jet final state (the two parton final state gives a contribution 
proportional to 5(0 — vr)) as 



Di 



as(/i) 



R ^^N "s(/i) , 



2. 

,\ r 

a,(0) + 27r6ofii,(0) log 

for 7^ TT , 



^2 

, CM , 



27r 



:2.10) 



(2.11) 



so that the expression for Gij becomes 



+ 



27r 



27r 60 log 



where 




127r 



+ 0{al), (2.12) 



(2.13) 



/i is the renormalization scale, Ga 
is the number of active fiavours. 



iVc = 3 for an SU(3) gauge theory, Tp = 1/2 and n/ 



■^Of course, this holds if we neglect terms that are suppressed by powers of m. Terms suppressed 
by powers of m and enhanced by powers of log(i?cM/™) may still (and in fact are) present in dj. 
Furthermore, uncalculable effects of order AqcdZ-Ecm rnay still be present. 

''A version of the program described in was used which implements the dipole subtraction method 
for massive partons |17| to combine real and virtual corrections. 
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We observe that expanding the ratio in powers of as is essential to implement the 
cancellation of large logarithms in Gij order-by-order in the perturbative series. If we had 
not done so, the numerator and denominator of Gij would contain large logarithms at 
some fixed order, and their ratio would be misleading (i.e. it would be sensitive to higher 
order terms with higher powers of logarithms) . 

In Fig.Qwe have plotted Gij at leading and next-to-leading order in Og, for i, j = 1, 2, 3, 
at a renormalization scale /i = E^m = 91.2 GeV, using as(-EcM) = 0.12, with Uf = 5 active 
flavours. The NLO curves are in good agreement with the SLD preliminary data. 

Conversely, the SLD data can be used to extract the value of as. In fact, from the 
knowledge of the numerical values of the coefficients of the as and contributions in 
Eq. ()2.12|) . we can fit Gij{(f)) to the data, keeping ag as a free parameter. 

We have computed tables for the ag and ag coefficients of Eq. (j2.12|) . that can be used 
in a more detailed analysis. These tables (and the programs to generate them) can be 
obtained from the authors. 

2.1 Secondary 6-quark and the running coupling 

We should mention here that, in the calculations of Ref. the b quark is treated as a 
heavy particle, and thus does not contribute to the running of the coupling constant. Thus, 
in principle, one should use this calculation in conjunction with the 4-flavour coupling 
constant ag^''. Since we are dealing with decays, however, secondary 6-quark pairs 
may in fact be produced (from gluon splitting), and when this happens, they behave, in 
some sense, as light flavours, contributing terms that effectively modify the running of ag. 

In Fig. we show the Feynman graphs for secondary heavy-flavour production, giv- 
ing rise to four massive-quark flnal state. The graphs in Fig. |21 also appear with the two 
identical quarks and/or the two identical antiquarks exchanged. When taking the ampli- 
tude squared, this gives rise to two interference terms. Their contribution is known to be 
extremely small, and will be neglected in the following. 

In addition, virtual 6-quark loops arise only in the self-energy insertion of the emitted 
(on-shell) gluon. Their contribution is zero in the renormalization scheme we use, since, 
in this scheme, the heavy flavour decouples for small momenta (see Appendix|B]for more 
details). Thus, the only effect of the heavy flavour is given by the graphs of Fig. |21 

In the quantities we are considering, the possibility of detecting secondary (instead 
of primary) B mesons has little impact, since secondary i?'s are very soft. If we make 
the approximation of neglecting effects that are suppressed by powers of the mass in the 
secondary b line, the contribution of the graphs of Fig. El (without interference) is easily 
included using the following prescription: 

- compute the cross section with the number of light flavours nif equal to 5; 

- use the 5 flavour coupling ag^"*; 

- do not include any other graph with secondary heavy flavours. 
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Figure 1 : Leading order (dashed) and next-to- leading order (solid) curves for Gij , 
i,j = 1,2,3, at a renormalization scale /U = E'cm = 91.2 GeV, using as{EcM) = 
0.12. The error bars in the preliminary data from the SLD Collaboration ^S] are 
the quadratic sum of systematic and statistical errors. 
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Figure 2: Secondary 6-quark production. 



In essence, the prescription amounts to including the secondary heavy flavour by simply 
increasing riif by one. When doing this, one in fact is also changing the renormalization 
scheme for the heavy flavour, by treating it as if it were light. Thus, the correct coupling 
constant to use in this framework is a'i'\ In Appendix^ we discuss this change of scheme 
in more detail and from a slightly different point of view. 

Secondary 6-quark production may also arise from primary light quarks, connected 
with the Z^/'-f vertex. We neglect altogether this mechanism, since it generates only soft 
b quarks. 



3 The m ^ limit 

As explained in Sec. El Gij{(j)) are free from potentially dangerous logarithms of the 
ratio Ecyi/m. In order to investigate the stability of our expressions in the m— >-0 limit, 
we introduce the integral of Gij over and we deflne (we drop the B /h subscript for ease 
of notation) 

= J ^ dcos0 = ^^.-i^^^i-i^ ' (3-1) 

where 

{^') = (4"') = A, (3.2) 



x'f'^xi'^) = I Ai(0)t^cos0. (3.3) 



-1 



To evaluate Eq. ()3.H) up to order a| and for i,j = 2, 3, we rewrite as follows 

{{l-x\-'){l-xi-'))-{l-x{-'){l-xi-') 
{x\ '){xi ) 
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and we expand in ctg 



{{I - X\-'){1 - Xi-')) : 



as(/i) 



' 2n 
Cij + 2tt bo bij log 



(3.5) 



^2 



2ti 



+ a 



(3.6) 



to obtain 



1 + 5..^ + 



Cjj + (oj + Oj) bij — ttittj + 2tx bo bij log 



^2 



CM / _ 



'tts(/^~ 

277 



(3.7) 

Note that, in Eq. ()3.4|) . the expectation value of ^ and X2 should be known only up 
to order ttg, and that, in the correlation term ((1 — x\~^){l — X2~^)), the two-jet final 
state gives zero contribution, since it is proportional to 5(1 — xi)6{l — X2), at all orders 
in as- This is the reason why we do not need to know the q;| two-loop virtual term for 
e+e-^QQ. 

The coefficients and bij can be computed analytically. We have collected their 
expressions in Appendix El In Tabled we give the numerical values for the coefficients 



Table 1: Results for the coefficients a^, bij and Cij appearing in the definition of Vij (see 
Eqs. fl3.5|) - ()3.7p ). The mass m is in GeV and /i = E^m = 91.2 GeV. 





m = 0.01 


m = 0.1 


m = 1 


m = 3 


m = 5 


m = 10 




29.1611 


20.9741 


12.7866 


8.87675 


7.05412 


4.56726 


as 


46.6475 


33.8553 


21.0574 


14.9209 


12.0317 


8.01007 




0.666666 


0.666589 


0.661895 


0.636326 


0.598452 


0.478637 


&32 


1.13333 


1.13321 


1.12569 


1.08448 


1.02307 


0.826778 


hs 


1.91111 


1.91093 


1.89945 


1.83563 


1.73915 


1.42415 


C22 


820.1(3) 


420.51(4) 


155.14(1) 


75.787(7) 


49.176(3) 


22.943(1) 


C32 


1290.5(4) 


663.99(7) 


247.19(2) 


122.15(1) 


80.053(6) 


38.268(3) 


C33 


2027.5(5) 


1046.5(1) 


393.09(3) 


196.60(2) 


130.23(1) 


63.871(5) 



Oj, bij and q^, for a down-type quark with mass m = 0.01, 0.1, 1, 3, 5 and 10 GeV, 
at a center-of-mass energy E^y^ = /i = 91.2 GeV. Observe that, while the individual 
coefficients and Cij can be quite large, due to the presence of logarithmic terms, the 
ratio Tij is well behaved, as illustrated in Table |21 This is due to the fact that the large 
logarithms present in the single coefficients cancel out in the ratio rij. 
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Table 2: Results for Vij. The mass m is in GeV and n = E^^ = 91.2 GeV. 



m 


^22 - 1 


r-i2 - 1 


r^-i - 1 


0.01 
0.1 

1 
3 
5 

10 


0.106103 as+0.215(4) ai, 
0.106091 as+0.217(l)af 
0.105344 as+0.2175(4) al 
0.101274as+0.2100(l)a2 

0.0952466 as+0.19905(9) al 
0.0761774 as+0.16351(4) al 


0.180375 as+0.403(6)a^ 
0.180356 as+0.407(2)a2 

0.179159 as+0.4066(7) al 
0.172600 as+0.3930(3) al 
0.162827as+0.3725(2) a| 
0.131586 as+0.30605(7) al 


0.304162 as+0. 750(9) ai, 

0.304133 as+0. 755(3) al 

0.302307as+0.7524(l)a2 

0.292149 as+0.7287(5) al 

0.276794^3+0.6919(3) al 
0.226661 as+0.5706(l)a2 



4 Theoretical uncertainties 

In this section, we investigate the theoretical uncertainties related to the renormaliza- 
tion scale and mass scheme used in the calculation. 

The renormalization-scale dependence of r is illustrated in Table El for the case m = 

5 GeV, using the two-loop evolution of as with as = 0.12 for fj, = Ecm = 91.2 GeV. For 
Ecyi/2 < /i < 2i?cM the variation of (r^j — 1) around the central value /i = £^cm is roughly 
±10% at leading order. It is reduced to ±5% at NLO. 



Table 3: Results for for different renormalization scales with as{Ec^) = 0.12, E^^ = 
91.2 GeV and m = 5 GeV. 





= 

LO 


-E'cm/2 
NLO 


= 

LO 


NLO 


= 

LO 


2 -Ecm 

NLO 


^22 - 1 
^32 - 1 
^33 - 1 


0.0127974 
0.0218776 
0.0371903 


0.014937(2) 
0.026116(3) 
0.045455(5) 


0.0114296 
0.0195393 
0.0332153 


0.014296(1) 
0.024903(2) 
0.043179(4) 


0.0103320 
0.0176629 
0.0300257 


0.013622(1) 
0.023667(2) 
0.040922(3) 



In order to investigate the mass-scheme dependence, we now explicitly distinguish 
between the pole mass m and the MS mass m(/i) (also called running mass). Using the 
relation between m and m(/i) 



m = m[jjL) < 1 H Lp 

27r 



2- -log 

2 ^ 



7^ 



+ O {at) 



(4.1) 



and defining 



m 



E, 



CM 



E, 



(4.2) 



CM 
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we have 



r-ij (z) = Tij {z (/i) + /\z (/i)) = rij{z{^)) + 



2ti 



Ah, {z{^)) + 



a. 



(4.3) 



with 



In 



dbij 
dz 



(4.4) 



We have evaluated Eq. ()4.3p for tTi^Ecm) = 3 GeV using the analytic expressions for the 
coefficients bij collected in Appendix |X| The value m^Ecyi) = 3 GeV corresponds roughly 
to what one gets by converting m at the 6-mass scale to m, and then evolving m up to 
/i = EcM- We have further studied the dependence of Eq. ()4.3j) on the renormalization 
scale /i. We have used the 2-loop evolution of m{fi), which gives m{Ec^/2) = 3.20 GeV 
and fn{2Ea^) = 2.83 GeV. The results for (fj^ — 1) are listed in Table El Using the 
MS definition for the mass, the theoretical uncertainties due to the scale ambiguity are 
practically of the same size as in the pole-mass scheme. 



Table 4: Results for Vij for different renormalization scales with as{Ec^) = 0.12, i^'cM = 
91.2 GeV and m{EcM) = 3 GeV. 





= 

LO 


-E'cm/2 
NLO 


= 

LO 


' EcM 

NLO 


= 

LO 


NLO 


^22 - 1 
^32 - 1 
^33 - 1 


0.0135360 
0.0230755 
0.0390735 


0.015444(2) 
0.026986(4) 
0.046904(8) 


0.0121529 
0.0207120 
0.0350579 


0.014884(2) 
0.025897(4) 
0.044811(7) 


0.0110323 
0.0187981 
0.0318083 


0.014265(2) 
0.024743(3) 
0.042675(5) 



Comparing the results for obtained with an on-shell mass m = 5 GeV with those 
obtained with an MS mass m{Ecyi) = 3 GeV, we found that, at leading order, the differ- 
ences between (r^j — 1) using the two different mass definitions amount to 5 -i- 6%. This 
theoretical uncertainty is reduced to less than 4% at NLO. 

An interesting question can be raised about the use of the pole mass versus the MS 
mass, that is to say, which one is more natural in this context. In practice, this is a very 
difficult question to answer. In inclusive processes, in general, the MS mass is better, since 
it accounts for large powers of logarithms of m/EcM multiplying powers of the mass at all 
orders in perturbation theory. The case we are considering, however, is not an inclusive 
process, since we are weighting the cross section with final-state kinematic variables. We 
are thus unable to give a full answer to this question. However, we can certainly ask 
whether, at the order we are considering, there is some numerical indication that using 
the MS mass accounts at least for a large part of the logarithmic terms in our expression. 
The mass dependence of the coefficients of rij is reported in Table El Calling r^'^^ the 
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K : Pole mass 

□ : runninig; mass 
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Figure 3: The quantities Xij, defined in Eq. fl4.5|l . plotted as a function of m, both in 
the pole-mass scheme and in the MS-mass scheme, with Ecm = = 91.2 GeV. 



coefficient of al in r^j, we consider the expression 

. ,4.5) 
as a function of logz. The same expression in the MS scheme is given by 



rjf(^)-rjf(0) , 4 



3 

1 - -log^ 



dh 



dz 



(z) , (4.6) 



where we use /i = -Ecm, and z = z{fj,). If the slope of the second expression versus logm 
is smaller than the slope of the first one, we can take this as an indication of the fact 
that it is more appropriate to use the MS mass instead of the pole mass. The results are 
displayed in Fig. El The value r^'^\o) has been taken to correspond to m = 0.1 GeV in 
Table 121 (the point m = 0.01 GeV has been excluded, because of the large errors). The 
points at m = 1, 3, 5 and 10 GeV are then plotted. It is quite clear from the figure that 
the logarithmic slope is smaller for the pole-mass scheme. Furthermore, the radiative 
correction to the mass term is always smaller for the pole-mass scheme than for the MS 
scheme. We thus conclude that the quantities rij seem to be more naturally described in 
the pole-mass scheme rather than in the MS-mass scheme. 
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5 Conclusions 



In this paper we have investigated new measurable quantities (related to the quantities 
introduced in Ref. 0) defined in e+e~ annihilation, that carry information about the 
angle-dependent B-B energy correlation. We have shown that these variables have a well- 
behaved expansion in terms of the strong coupling constant a^, since large logarithms in 
the ratio of the mass of the final quarks over the centre-of-mass-energy cancel out order 
by order, together with non-perturbative effects of order Aqcd/^- 

We have compared the computed NLO results with the available preliminary data from 
the SLD Collaboration ^j, and found good agreement. In addition, we have generated 
tables ^ that can be used in a fit to the data, in order to extract the strong coupling 
constant at next-to-leading order. 

We have investigated the renormalization-scale and mass-scheme dependence of our 
results, and found that, varying the renormalization scale fi in the range -Ecm/2 < fi < 
2EcM, the scale dependence of (r^j — 1) is reduced from roughly ±10% at leading order to 
±5% at next-to-leading order. The theoretical uncertainties related to the mass definition 
are less than 4% at NLO. We have shown that these quantities seem to be more naturally 
described in the pole-mass rather than in the MS-mass scheme. 

From the discussion given in this paper, one may wonder if similar quantities may be 
defined for massless quarks, given the fact that mass singularities cancel in the variables 
we are considering, and thus they remain well defined in the massless limit. This is not 
quite the case. In fact, it is crucial for our discussion that the secondary heavy-fiavour 
production may be neglected. If this were not the case, evolution of the fragmentation 
function would become multi-dimensional, involving several fiavour components and a 
singlet contribution. Thus, the advantage of having a heavy flavour in the final state lies 
in the fact that the evolution becomes essentially non-singlet, so that we are left with 
a single component of the fragmentation function, that simply cancels in the ratios we 
propose. 
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A Analytic results for a2 5 0^35 ^225 ^325 ^33 

In this appendix, we collect the analytic results for the coefficients and bij, i,j = 2, 3, 
defined in Eqs. ()3.5|) and ()3.6|) . Calling the generic coefficient F{z), where z = m? / s and 

''The programs and the tables can be obtained upon request from the authors. 
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s — E'cM) we write 



2C 



where /3 = y/T^^4z and a; = (l— + The electro weak couphngs cy and are 
given by 

cv = Qlr + 2g'yQ,Re[x{s)]r' + g'y'\x{s)\'f', 

CA = y^'|x(^)|V^^ (A.2) 

with 

r = i-A_A+, 

/^^ = (l-A_A+)((?^2 + ^?f)-2(A_-A+)5e^5^, 

= -(l-A_A+)(7^ + (A_-A+)^7l. (A.3) 

Here, Q5 = —1/3 is the electric charge of the bottom quark, and g^/y denote the axial 
and vector couplings of the fermion /. In particular, gy = — | + 2sin^'(9vi/, S'a = ^| for an 
electron, and gy = + ^sin^i^^w, g\ — ~\ a bottom quark, where is the weak 
mixing angle. The function reads 

^^^^ ^ 4sin2i?vFCOs2i?vF s-Ml^iMzYz' ^^'^^ 

where Mz and Y z stand for the mass and the width of the Z boson. Finally, A_ (A+) 
denotes the longitudinal polarization of the electron (positron) beam. 

A simple calculation using the Born differential cross section for e^e~ -^bbg gives 

a^{z) = (-l + z-^zAln{u;)+(-^ + lz + 7zA(^, (A.5) 



6 

' ' ^ (^--10^+12^2^-^=^jln(a;) + (^--y^-10^2j|, (A.6) 



z 



25 10 5 , 32 11 

al[z) = [-— + —z + -z^-—z^-—z'']ln(uj) 



32 3 11 4A w 

— z z In t 

24 3 2 3 12 y ^ 



/ 433 137 779 2 3\ ^ , , 

+ \ ^^ + 112;^ —, (A.7) 

V 12 2 6 y 24' ^ ^ 

ai{z)-a^{z) _ /55 , , , 40 , , 15^, 



^2 - 20^ + 21^^ + —z^ + —z^j ln{uj) 
36 18 6 7 4' ^ ^ 



bUz) = 2z{l + z')Hu;) + (^^ + ^z-^z'-z'^f3, (A.9) 
fe^2(^) - bUz) ^ 2z{-2 + 3z-5z')ln{u;) + (^-l-^Jz + lz' + 5z')p,{A.lQ) 
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.17 1721 47 2 103 3 1 A /3 

+ \ z^ z^ z^--z^]-, (A.ll) 

\20 180 45 18 3 7 2 ^ ^ 



bi,{z)-b^, 



= ;2(-7+10^-10^=^-14/)lnH 

/ 27 287 16 2 37 o „ A ^ 

+ z + —z^ + —z^ + 7z^] A.12 

V 20 60 15 6 r ' ^ ' 

I.V , ^ /14 13 „ 2 32 3 16 4 5\, , , 

^Uz) = z{^—-—z-Sz^ + —z'' + —z^-z''yn{u) 

43 , 533 221 ^ 1037 3 31 4 5^ /3 



— + —Z - — - — + -, (A.13) 

30 45 36 90 6 / 2' ^ ' 

hUz) - bUz) _ ^ ^ ^2 + 17^ + -z^ - I2z' - 1%z^ - -z^) ln(o;) 

3 3 / 

109 661 , 53 2 157 3 287 4 A ^ 



+ 2 + — 2^ + 2^ + 2^ + 35^^-. (A.14) 

V 15 30 20 6 6 / 3 ^ ^ 

In the m — i> hmit, we have /? — > 1, ^ 2; = vn} j s — 0, so that 

«2 = 2C^(-^ln(^)-H) , 

as = 2Cp(--\m(z)- — 
^ V 24 ^ ^ 288 

622 = 2Cf- , 

633 = ICf^. (A.15) 
60 

This behavior is well illustrated from the collections of values in Table ^ while the bij 
are finite in the m— i>0 limit, the coefficients diverge as log(m/i?cM)- 



B Renormalization schemes 



Our calculation was carried out in the "mixed" renormalization scheme of Ref. 
in which the light flavours riij are subtracted in the MS scheme, while the heavy-flavour 
loops are subtracted at zero momentum. In this scheme the heavy flavour decouples at 
low energy. In fact, convergent heavy-flavour loops are suppressed by powers of the mass 
of the heavy flavour. The only unsuppressed contributions come from divergent graphs. 
But those are subtracted at zero external momenta, so their contribution is removed by 
renormalization for small momenta. In the mixed scheme, charge renormalization is given 
by the prescription 



1 - as(yu)T 




+ O (at) 



(B.l) 
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where al is the bare couphng constant, and as{fi) is the mixed-scheme couphng constant 
at the scale fi. We have defined 

127r 

and 

- = - + log(47r)-7s, (B.3) 
e e 

where dimensional regularization is used, with the number of space-time dimensions equal 

to 4 - 2e. 

In the MS scheme the renormalization prescription is 



3^ = /i2^a3(^)(l-«s(^)-&^'^| , (B.4) 



where nj = riif + 1, and as(/i) is the MS coupling constant at the scale fi. The cou- 
plings in the two schemes obey a 4-fiavour and 5-fiavour renormalization-group equation 
respectively 

r«s(^) = -bo''^'^ alifi) , ^ n ftsM = -b'-Q^^alifi) , (B.5) 



c/log/i^ ^ ' dlogyU^ 

that can be easily derived by imposing the constancy of the bare coupling under a 
renormalization group transformation. Combining Eqs. (jB.lll and ()B.4j) we have 

as(/i) = «s(/i) - ^ log (^^^ alifi) + O (af) , (B.6) 

which is the standard MS matching condition for fiavour crossing. 

In the heavy fiavour production calculation we are considering, if we express the result 
(usually given in terms of a^) in terms of ctg, an extra term arises, equal to 



Bornx«s(/x)^logf4) • (^-7) 



The only other diagrams where secondary heavy fiavours appear are the graphs with a 
gluon splitting into a heavy-fiavour pair. It can be shown that, if we neglect powers of 
m/EcM, the term in Eq. ()B.7|) is exactly the term one needs to turn the calculation of 
the heavy-fiavour splitting, which is regulated in the collinear limit by the heavy-fiavour 
mass, into the corresponding MS subtracted calculation, provided that interference terms 
of the heavy fiavours coming from gluon splitting with the primary heavy fiavours are 
fully neglected. 

Equation ()B.7|) follows immediately from the theory of heavy-fiavour fragmentation 
functions. According to Ref. ^^1; Eq- (3.16), the cross section for the inclusive production 
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of a heavy quark via gluon splitting, with a fraction x of the energy of the gluon, is (at 
leading order and neglecting powers of tti/Ecm) 



da ctsTp 
dx ° 27r 



x^ + (l-xf 



log 


(4) 







da , r, , 



(B.8) 



where cxo is the cross section for the production of the gluon and a stands for the massless 
limit, MS subtracted cross section for the gluon splitting process. Integrating both sides 
in X one finds 

asTp , / /i^ \ 
a = a- (To— — log \^] ■ 



(B.9) 



The second term on the right-hand-side is precisely the contribution of Eq. ()B.7|) . 
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